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ABSTRACT 

We present Space Telescope Imaging Spectrograph observations of 14 nearby low-luminosity active galactic 
nuclei, including 13 LINERs and 1 Seyfert, taken at multiple parallel slit positions centered on the galaxy 
nuclei and covering the Ha spectral region. For each galaxy, we measure the emission-line velocities, line 
widths, and strengths, to map out the inner narrow-line region structure, typically within ^100 pc from the 
galaxy nucleus. There is a wide diversity among the velocity fields: in a few galaxies the gas is clearly in 
disk-like rotation, while in other galaxies the gas kinematics appear chaotic or are dominated by radial flows 
with multiple velocity components. In most objects, the emission-line surface brightness distribution is very 
centrally peaked. The [S II] line ratio indicates a radial stratification in gas density, with a sharp increase within 
the inner 10-20 pc, in the majority of the Type 1 (broad-lined) objects. The electron-density gradients of 
the Type 1 objects exhibit a similar shape that is well fit by a power law of the form n e = «o(r/l pc)", where 
a = -0.60 ±0.1 3. We examine how the [N II] A6583 line width varies as a function of aperture size over a range 
of spatial scales, extending from scales comparable to the black hole's sphere of influence to scales dominated 
by the host galaxy's bulge. For most galaxies in the sample, we find that the emission-line velocity dispersion 
is largest within the black hole's gravitational sphere of influence, and decreases with increasing aperture size 
toward values similar to the bulge stellar velocity dispersion measured within ground-based apertures. We 
construct models of gas disks in circular rotation and show that this behavior can be consistent with virial 
motion, although for some combinations of disk parameters we show that the line width can increase as a 
function of aperture size, as observed in NGC 3245. Future dynamical modeling in order to determine black 
hole masses for a few objects in this sample may be worthwhile, although disorganized motion will limit the 
accuracy of the mass measurements. 

Subject headings: galaxies: active - galaxies: bulges - galaxies: kinematics and dynamics - galaxies: nuclei 



1. INTRODUCTION 

With its ~ 0"1 resolution, the Space Telescope Imag- 
ing Spectrograph (STIS) aboard the Hubble Space Telescope 
(HST) has provided the exceptional ability to spatially resolve 
gas kinematics within the gravitational sphere of influence 
of the putative central black hole in nearby galaxies. Re- 
solving the gas kinematics very close to the black hole al- 
lows for the determination of black hole mass in a relatively 
straightforward manner, provided that the gas is in Keple- 
rian rotation in a thin, disk-like structure. STIS has thus 
had a major impact on the field of supermassive black hole 
detection in galactic nuclei, and data from STIS have con- 
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tributed greatly toward establishing the correlations between 
black hole mass and host-galaxy properties, such as those 
with the stellar veloci t y dispersion ( Ferrarese & Merritt 2000; 
Gebh ardt et al.l 120001 iTremaine et alj|2002l) and bulge lumi- 
nositv dKormendv & Gebhardtll2001l) . 

More recently, ground-based integral field unit (IFU) obser- 
vations have provided two-dimensional (2D) kinematic infor- 
mation on the ioniz ed gas in galaxy centers (Sarz i et aT1l2006l : 
iDumas et al.ll2.Q07b . Additionally, ground-based IFU observa- 
tions with the assistance of adaptive optics h ave produced gas 
dynamical black hole mass measurements dHicks & M alkan 
120071: fNeumaver et al. 2007). Without the use of adaptive op- 
tics, however, ground-based observations are usually unable 
to resolve scales comparable to the black hole radius of in- 
fluence in nearby galaxies. Although STIS is a long-slit in- 
strument, 2D information similar to IFU data can be obtained 
with STIS by using multiple slit positions. 

Even with the past decade of high-resolution observations, 
only about 30 dynamical black hole mass measurements have 
been made (Fer rarese & Fordl2005l) . Consequently, searching 
for regular gas velocity fields dominated by gravitational mo- 
tion is essential for extending the number of black hole mass 
measurements and further establishing the connections be- 
tween black holes and host galaxies. While searching for reg- 
ular gas velocity fields is clearly important, it appears, how- 
ever, that they are fairly rare. Instead, complicated and chaotic 
emission-line velocity fields are more commonly found nea r 
the centers of galaxies dHo et alJl200llAtkinson et al.ll2005[) . 
Although only upper limits on the black hole mass can be 
estimated using gas dynamical modeling in these situations 
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dSarzi et al.ll2002f) . subarcsecond-resolution observations can 
still be useful for studying the dynamical state of the gas 
and searching for radial motions that might be related to ac- 
tive galactic nucleus (AGN) fueling or outflows. Modeling 
AGN outflows on subarcse cond scales (ICapetti et al.l Il996t 
ICrenshaw & Kra emer 2000) is important for furthering our 
understanding of the connection between AGN feedback and 
host-galaxy properties. Of particular interest are the connec- 
tion to star formation, the growth of the supermassive black 
hole, and how gas may be driven out f rom the host galaxy 
as predicted in some feedback scenarios dSpringel e t al. 2005; 
iHopkins et al1l2lM[Sr^cMetini2007h . 

Moreover, the narrow-line region (NLR) velocity disper- 
sion ((T g ), often measured from the width of the [O III] A5007 
line, may provide a reasonable estimate of the host-galaxy 
stellar velocity dispersion (a*), albeit with su bstantial scatter 
dNelson & Whittled 1996tlGreene & Holl200l . Detailed NLR 
studies are essential in understanding the origin of the scatter 
in the er^-ov relationship. This relationship can be useful for 
distant and luminous AGNs where mea surement of the stel- 
lar velocity dispers ion is not feasible dSalviander et al.ll2007l; 
IShields et al.ll2003l Using the narrow emission-line widths 
as a proxy for stellar velocity dispersion is only meaningful if 
the NLR gas is dominated by virial motion in the host-galaxy 
bulge and not by the gravitational influence of the black hole 
or nongravitational motions. Therefore, outflows and similar 
nongravitational motions, which can disrupt viral motion, are 
one potential cause of the scatter in the Cg-cr* relationship. 

While nongravitational motions are clearly important in 
setting the narrow-line widths in AGNs with strong out- 
flows, there can be an important contribution to the line 
widths from nongravitational motion even in objects that 
are not outflow-dominated . From HST STIS observations, 
IVerdoes Kleijn et all d2006l) found that circumnuclear gas in 
Fanaroff-Riley 1 radio galaxies typically has an intrinsic ve- 
locity dispersion in excess of that expected from purely grav- 
itational motion in the host-galaxy potential. They concluded 
that this excess line width is driven by the injection of energy 
by the AGN, although there was no obvious relationship be- 
tween the AGN radiative luminosity and the magnitude of the 
nongravitational line dispersion. 

Oth er origins for the scatter were studied by iGreene & Hoi 
(2005) using a large sample of Sloan Digital Sky Survey 
(SDSS) type 2 AGNs. They looked for correlations of nuclear 
properties (radio power, AGN luminosity, Eddington ratio), as 
well as global properties (host-galaxy morphology, local envi- 
ronment, star-formation rate), with the excess line width, de- 
fined to be the difference between the dispersion of the [O III] 
A5007 line and the stellar velocity dispersion. They found a 
strong correlation between the excess line width and Edding- 
ton ratio, where larger excess line widths, indicative of out- 
flowing gas, were seen in the higher Eddington ratio objects. 

Another potential explanation for the scatter in the a g - 
relationship i s the effect of observational aperture size. 
iRice et all d2006l) investigated the relationship between line 
width and aperture size in a sample of nearby Seyferts us- 
ing single-slit STIS observations. They found a wide range 
of behavior in how the widths of the [O III] A5007 and [S II] 
A A67 1 6, 673 1 lines varied as a function of radius, but the most 
common trends were either a roughly constant width or an in- 
crease in width with increasing aperture size. Their analysis 
showed that the line widths measured within the largest STIS 
aperture were systematically smaller by 10%-20% than both 
the stellar velocity dispersion and the line widths measured 



within a ground-based sized aperture. They concluded that 
the smaller line width measured within the randomly oriented 
STIS slit is the result of sampling a smaller portion of the ve- 
locity field than the larger ground-based aperture. 

Beyond the cr g -<J* relationship, we can addre ss oth e r open 
questions concerning low-luminosity AGNs. lLaorl (2003) 
pointed out that while the NLR in luminous AGNs occurs on 
scales of tens of parsecs, where its dynamics are dominated by 
the bulge, the NLR in low-luminosity AGNs is more compact, 
and its dynamics could be dominated by the black hole. It is 
therefore interesting to ask whether the narrow-line widths in 
low-luminosity AGNs are set by the black hole or by the bulge 
potential. 

The NLR emission lines also provide information on 
ionization mechanisms. Spectr a of low-ionizat ion nuclear 
emission-line regions (LINERs; lHeckmanlfT980h can be re- 
produced with photoionization models invol ving a nonstellar 
continuum and low ionization parameters dFerland & Netzerl 
I1983L lHalpern & Steined 119831) . but high electron densi- 
ties, as well as a large range of densities, are needed 
in o rder to reproduce the strength of the \Q III] A4363 
line dFilippenko & Halperrill984]Fflippenkd[l985l: iHo et ail 
119931) . iFilip penko & Halpern ( 1984J) demonstrated that such 
conditions exist in the LINER NGC 7213, and that there is 
a correlation between line width and the critical density for 
collisional deexcitation. A similar relationship between line 
width an d ionization pote ntial was also previously seen in 
Seyferts dPelatet al.lll981l) . This implied that the NLR is not 
a homogeneous environment, but must have radial gradients 
in velocity and density. In ground-based spectra these spatial 
gradients were generally unresolved, but using STIS it is pos- 
sible to directly resolve the density grad ients in the inner NLR 
dBarth et al.ll2001at IShields et al.ll2007l) . It is therefore of in- 
terest to determine what fraction of LINERs show spatially 
resolved density gradients, and whether the type 2 LINERs, 
where a central photoionizing source is not always seen, ex- 
hibit density gradients. 

In this paper, we present multislit HST STIS observations 
of 14 low-luminosity AGNs, including 13 LINERs and 1 
Seyfert. We use the high spatial resolution of STIS to map 
out the 2D kinematic structure of the NLR within the inner 
~ 100 pc, and search for velocity fields that show signs of 
regular rotation for which gas dynamical modeling may be 
performed. The subarcsecond resolution STIS data also al- 
low us to detect electron-density gradients through the [S II] 
A6716/A6731 line ratio diagnostic. Finally, we examine how 
the [N II] line width varies as a function of aperture size over 
spatial scales ranging from those comparable to the black hole 
sphere of influence to scales dominated by the bulge. 

2. OBSERVATIONS AND DATA REDUCTION 

Our goal is to examine the NLR of low-luminosity AGNs 
through STIS spectra of the Ha spectral region at multiple slit 
positions. Eight galaxies observed under program GO-7354 
(PI Dressel) fit this description, and we retrieved those from 
the HST archive. We combined these eight galaxies with six 
galaxies that were observed as part of our own program, GO- 
7403, to arrive at a sample of 14 galaxies. Program GO-7354 
was originally designed to study LINER nuclei with compact, 
flat-spectrum radio cores, while the observations under pro- 
gram GO-7403 focused on obtaining black hole masses in 
a sample of low-luminosity broad-lined AGNs. All of the 
galaxies are optically classified as LINERs except NGC 3227, 
which contains a Seyfert 1 nucleus, and NGC 3245, which is 
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a transition object. Properties of all 14 galaxies in the sam- 
ple are given in Table Q] and details of the observations are 
presented in Table [2] 

With the exception of NGC 3998, NGC 4594, and NGC 
6500, all of the galaxies from program GO-7354 were ob- 
served with the STIS 52x0.1 aperture with no gap between 
the three adjacent slit positions. NGC 3998, NGC 4594, and 
NGC 6500 were observed with the same 52x0 . 1 aperture, 
but consist of five slit positions, again with no gap between 
neighboring slit positions. The CCD was read out with bin- 
ning by two along the dispersion axis, creating a scale of 
0"0507 pixel" 1 along the spatial axis and 1.108 A pixel" 1 
along the dispersion axis. The slit position angle was con- 
strained to be within 15° of the major axis of the gaseous 
disk. The G750M grating provided coverage of 6490-7050 
A, which included the [N II] AA6548, 6583, Ha, and [S II] 
AA6716, 6731 lines. For this program, a typical exposure time 
per position was ~ 400 s. 

Observations of the six galaxies under program GO-7403 
were obtained with the STIS 52x0.2 and a spacing of 0."25 
between adjacent slits (i.e., a gap of 0."05 between slits). Two 
of these galaxies, NGC 1052 and NGC 3227, were observed 
with seven parallel slit positions, while NGC 3245, NGC 
4036, NGC 4278, and NGC 4579 were observed with five 
slit positions. For these galaxies, the CCD was read out in 
unbinned mode, and the slit position angle was left uncon- 
strained. With this setting the 2D spectral image has a scale 
of 0"0507 pixel" 1 along the slit and 0.554 A pixel" 1 along the 
dispersion direction. The G750M grating was used to cover 
6300-6860 A, which included the Ha+[N II], [S II] AA6716, 
6731, and [O I] AA6300, 6364 emission lines. A typical ex- 
posure time per position was ~ 3000 s. The combination of 
longer exposure times and the use of a wider slit made it pos- 
sible to map out the NLR to larger radii and at higher signal- 
to-noise ratio (S/N) than for program GO-7354. 

The data for all 14 galaxies were reduced using the standard 
STScI pipeline. Dark and bias subtraction, as well as flat- 
field corrections, were performed, along with the rejection of 
cosmic -ray events through the combination of multiple subex- 
posures. The data were wavelength and flux calibrated and 
rectified for geometric distortions. Occasionally a number of 
hot pixels remained even after the data were run through the 
pipeline. In those cases, we performed an additional clean- 
ing step by replacing the bad pixel with a local mean value in 
the flattened data, and then applied the geometric distortion 
correction to the cleaned images. 

In addition to the STIS observations, we obtained Wide 
Field Planetary Camera 2 (WFPC2) images from the archives 
in order to examine the nuclear morphology. When available, 
we retrieved the images from the WFPC2 Associations web 
site 9 . Otherwise, they were retrieved from the HST archive. 
In order to remove cosmic rays, "CR-split" observations were 
combined using the IRAF combine task. If only single 
exposures were available, they we re cleaned using the LA- 
COSMIC task (Ivan Dokkumll200Tl) . 

Black hole mass measurements have been previ ously made 
for thre e of the galax ies in the sample: NGC 3 245 (Barth et al. 
2001b ). NGC 3998 (Ide Francesco et al.l2006h. andNGC 5077 
dde Francesco et al. 1 120081) . In addition, de Fra ncesco et al.l 
(2008) presented velocity curves for NGC 6500 and IC 989; 
however, both velocity fields were too complex to model the 

9 http://archive.stsci.edu/hst/wfpc2/search.html . 



gas as a thin circularly rotating disk, so no black hole masses 
were determined for these two galaxies. 

3. MEASUREMENT OF EMISSION LINES 

Spectra were extracted from individual rows of the geomet- 
rically rectified 2D images, and the rows were extracted as far 
out as emission lines were detectable. For the eight galaxies 
from program GO-7354, the emission lines were typically de- 
tectable out to about 0"5 from the slit center, and for the six 
galaxies from program GO-7403, emission lines were mea- 
sured out to about l."5. Before any fitting was performed, the 
continuum was removed from the spectrum on a row-by-row 
basis by fitting a straight line to the continuum regions cov- 
ering 6280-6380 A, 6500-6650 A, and 6700-6750 A. This 
continuum fit was then subtracted from the spectrum. In gen- 
eral, a straight line is not a very good continuum model near 
Ha, but given the low S/N in the continuum at most positions, 
the small wavelength range, and the AGN dominance in the 
nuclei of some objects, a proper starlight subtraction could 
not be performed. During the emission-line fitting process, 
multiple rows that contained weak emission lines far from the 
slit center were binned together to improve the S/N. 

In order to obtain the emission-line velocity, velocity dis- 
persion, and flux, we fit each spectrum with a set of Gaus- 
sians, where the fit was carried out row-by-row with a 
Levenberg-Marquardt least-squares minimization technique 
using C. Markwardt's "MPFIT" library in IDL 10 . Single 
Gaussians were simultaneously fit to all the emission lines 
present in the spectra that were strong enough to be fit accu- 
rately. The fluxes of the [N II] AA6548, 6583 lines were held 
at a 3: 1 ratio, and all lines were required to have a common ve- 
locity. The [N II] lines were constrained to have equal velocity 
widths. The widths of the [S II] lines were also required to be 
the same, which was necessary due to the lower S/N in the 
off-nuclear rows. Although the widths of the two \S III lines 
can differ, the effect is sm all (IFilippenko & Halperni Il984t 
iFilippenko & Sargenj[l988l) . If present, the [O I] line widths 
were constrained to be equal. 

Whenever this basic model did not return adequate fits due 
to either low S/N or to severe blending between the [N II] 
and Ha lines, we applied further constraints, such as tying to- 
gether the Ha and [N II] widths, or linking the Ha, [N II], 
and [S II] widths. Additionally, some spectra appeared to 
contain two kinematically separate components at some lo- 
cations. In these instances, we fit two Gaussians to each of 
the Ha, [N II] AA6548, 6583, [S II] AA6716, 6731, and the 
[O I] A6300 lines. We fit a single Gaussian to the [O I] A6364 
line, and sometimes to the [O I] A6300 line, in place of two 
Gaussians since often these lines were too weak to constrain 
a multiple-component model. Each set of Gaussians in the 
multiple-velocity component model was constrained in the 
manner previously described above for the basic model. 

Near the nuclear regions, several galaxies in the sample ex- 
hibited a broad Ha component in addition to the narrow lines 
of Ha and [N II] AA6548, 6583. For these regions, an addi- 
tional component in the form of either a single or a double 
Gaussian was added to the fit. If the broad component could 
not be adequately fit by two Gaussians, we selected points 
along the broad base of the Ha line and fit a cubic spline 
to trace the shape of the broad component. This was then 
subtracted from the spectrum, leaving just the narrow emis- 
sion lines. A cubic-spline model was needed for several rows 

10 http://cow.physics.wisc.edu/~craigm/idl/idl.html . 
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while fitting the broad component in NGC 3227, NGC 4278, 
and NGC 4579. We judged the effectiveness of an additional 
Ha broad component both by eye and by comparing the re- 
duced x 2 , where if the reduced % 2 was improved by 15% or 
more, we accepted the broad-component fit. Figure Q] shows 
examples of fits with single or double narrow-line compo- 
nents and with a broad Ha component. 

From the model fits, we obtained the velocity, velocity dis- 
persion, and flux for the [N II] A6583 emission line. In cal- 
culating the velocity dispersion, we accounted for the instru- 
mental line width by using observations of the STIS wave- 
length calibration lamp. Gaussian models were fit to the 
emission-line profiles of extracted calibration-lamp spectra. 
The variation in line width as a function of wavelength was 
insignificant across the spectrum and as a result, we used the 
average value from measurements of several lines. The instru- 
mental velocity dispersion was er,- w 3 1 km s" 1 for program 
GO-7403 and a-, w 17 km s _1 for program GO-7354. These 
instrumental line widths were then subtracted in quadrature 
from the observed line widths. 

4. RESULTS 

The velocity, velocity dispersion, and flux measured from 
the [N II] A6583 emission line for 12 of the 14 galaxies is plot- 
ted as a function of the position along a single slit in Figures 
l2al42nl We show the galaxy's systemic velocity, taken from 
optical line measurements in the Third Reference Ca talogue 
of Bright Galaxies (RC3) dde Vaucouleurs et al.|[l991l) . as the 
dashed line in the radial-velocity plots. For each galaxy, we 
also calculated the flux-weighted average velocity measured 
from the largest square aperture possible given the number of 
STIS slits for the galaxy and the number of rows from which 
spectra could be extracted. For five galaxies in the sample, the 
flux-weighted mean velocity differs from the RC3 velocity by 
>60kms _I . 

We define the galaxy center (y-offset = 0) to be the row with 
the largest continuum flux. The peak of the continuum flux 
coincided with the peak of the Ha broad-component flux in 
five of the seven broad-line galaxies, and in the cases of NGC 
4579 and NGC 3998, the peak broad-component flux was off- 
set from the peak of the continuum by one CCD row, corre- 
sponding to 4 and 8 pc, respectively. One possible explanation 
for the small offset is differential extinction between the con- 
tinuum and the emission-line flux. Similarly, the peak of the 
narrow-line emission coincided with the peak of the contin- 
uum in five of the seven narrow-line galaxies, but for NGC 
4036 and IC 989, the peak of the [N II] A6583 flux was off- 
set from the continuum peak by one CCD row, corresponding 
to 5 and 50 pc, respectively. The Ha broad-component flux 
peak coincided with the narrow-line flux maximum in only 
two (NGC 1052 and NGC 4278) of the seven type 1 AGNs. In 
the other type 1 objects, the offset between broad and narrow- 
flux peaks was just one CCD row. 

For completeness we include the velocity fields from 
NGC 3245, NGC 3998, and NGC 5077, although th e 
data have been previously pre sented by iBarth et al.l d2 001b). 
Ide Francesco et al.l (12006!). and lde Francesco et al.l (120081) . re- 
spectively. We a lso include the velocity fie lds for NGC 6500 
and IC 989 since lde Francesco et al.l(l2008l) present the veloc- 
ity curves from the central slit positions, while we display the 
velocity fields from every slit position. Our measurements of 
the velocity fields from all five of these galaxies are consistent 
with those results previously published. 

For all but two galaxies, a continuum-subtracted Ha+[N II] 



image was created by subtracting a scaled WFPC2 702W im- 
age, or a scaled combination of F547M and F791W images, 
from a F658N image. This is shown along with the positions 
of the STIS slits in the second panel of the set of WFPC2 im- 
ages at the top of the figures. We verified that the central slit 
position was placed on the nucleus by examining the peak-up 
image, which is an undispersed image taken through the slit, 
and the acquisition image, which is an initial image of the 
object taken with a larger aperture. We compared the spatial 
profile of the peak-up image to the spatial profile produced by 
vertical cuts through the nucleus of the acquisition image, and 
found similar profiles, indicating that the slit was centered on 
the nucleus. 

The 14 galaxies in the sample show a large diversity in their 
observed velocity fields. Some velocity fields are consistent 
with rotation but perhaps at a low S/N, some galaxies seem 
to be in overall rotation with significant random components, 
some appear to be dominated by irregular motions, and others 
demonstrate clear evidence for outflows. As mentioned ear- 
lier, three galaxies (NGC 3245, NGC 3998, NGC 5077) have 
velocity fields previously shown to be consistent with circu- 
larly rotating disks. While there is no foolproof way to predict 
which galaxies will contain regularly rotating gas, iHo et al.l 
d2002l) determine that the morphology of dust lanes can be 
a reliable indicator of organized velocities fields. They find 
that galaxies with smooth, circularly symmetric dust rings of- 
ten have velocity fields suitable for gas dynamical modeling, 
while those objects with clearly disorderly dust lanes almost 
never contain regular velocity fields. 

Figure [3] shows the [N II] A6583 surface brightness as 
a function of projected angular distance from the nucleus, 
where each point corresponds to a single row measurement. 
In order to emphasize the trends, we averaged together points 
within 0."l-wide annuli for the galaxies observed under pro- 
gram GO-7403. For NGC 3998, NGC 4594, and NGC 6500, 
we averaged together points that fell within 0."05-wide bins. 
These averages are plotted as black triangles. The shape of 
the surface brightness profile varies between galaxies, where 
some galaxies show a dramatic increase in surface bright- 
ness within the innermost regions and other galaxies exhibit a 
gradual rise in surface brightness over a larger range of radii. 
However, almost all of the galaxies in the sample do exhibit a 
centrally peaked surface brightness profile. NGC 3998 shows 
the sharpest peak in surface brightness, where the [N II] sur- 
face brightness falls by about a factor of 58 between r = 
and r = 0."28. Conversely, NGC 4278, NGC 3245, and IC 
989 show the flattest profiles, where the [N II] surface bright- 
ness decreases by a factor of about 3-4 within the inner 0."25. 
NGC 2911 does not exhibit a central peak in the emission- 
line surface brightness, although the error bars are large and 
may mask a weak rise in surface brightness at smaller radii. 
In general, both type 1 and 2 LINERs have central peaks in 
emission-line surface brightness. 

Radial-velocity maps of each galaxy are shown in Figure[4] 
For galaxies having multiple velocity components, we show 
the radial-velocity map for the velocity component plotted 
in black in Figures [2a] |2dJ |2hl and |2l] In these multiple-slit 
velocity curves for NGC 1052, NGC 3227, NGC 4278, and 
NGC 4579, we attempt to label what appears to be coher- 
ent kinematic systems by using black and grey colors. How- 
ever, we cannot definitively say that all black points represent 
a single component and all grey points represent a separate 
component. This ambiguity in the identification of coherent 
velocity components affects the radial-velocity maps in Fig- 
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ure[4]and may cause the velocity field we show to be more 
chaotic than it is in reality. While the multiple-slit velocity 
curves presented in Figures I2aj - l2n1 contain the same informa- 
tion, the velocity maps show more clearly whether the gas is 
dominated by rotational motion or chaotic motions. 

Below we discuss the velocity fields for each galaxy in the 
sample, with the exception of NGC 3245 and NGC 3998. 

4.1. NGC 1052 

Overall the gas does not exhibit a regular rotation curve, 
and the presence of strong outflows can clearly be seen in 
the STIS data. The multiple velocity components identi- 
fied in our spectra are spatially coincident with the structure 
seen in the narrow -band Ha+[N II] emission-line image from 
Pogge ^t~ai1(l2000l ). who detected an ionization cone at a posi- 
tion angle of 96°. Two kinematically distinct components are 
seen in four of the off-nucleus slit positions and correspond 
to the location of jet-like features seen in the WFPC2 F658N 
image. The velocity splittings are as large as ~ 400-600 km 
s" 1 . The jet-like feature emerging from the nucleus toward 
the east is blueshifted by ~ 560 km s" 1 relative to the galaxy's 
systemic velocity of 1474 km s" 1 , and the western jet is red- 
shifted by ~ 240 km s" 1 . In the eastern jet, in slit position 7, 
the line widths are extremely large, with a velocity dispersion 
of 60010ns- 1 . 

Prev ious ground-based spectroscopic work by Sug ai et all 
(2005) has also found two components to the emission-line 
profiles, where one component is due to disk rotation at radii 
greater than ~ 1" with an axis of rotation at a position an- 
gle of 153°. The second component is attributed to a bipo- 
lar outflow, seen at a radius of ~1" from the nucleus, whose 
axis is at a position angle of ~ 103°. The blueshifted side 
of the outflow has a higher velocity of at least -1200 km 
s _1 relative to the systemic velocity, while the redshifted side 
has a velocity of 800 km s _1 relative t o the systemic veloc- 
ity (ISugai et alJ l2005|). Additional ly, iDavies & Illingworthl 
(119861) and lSchulz & Henkell (120031) found evidence for out- 
flows on large scales from ground-based spectra. 

The fit to the spectrum from the nuclear region of the central 
slit was improved by the addition of a broad Ha component, 
as expected since this is a well-known type 1.9 LINER; the 
broad Ha emission has been shown to be polarize d, indicating 
that th e broad-line region is seen in reflected light (Bart h et all 
fl999h . 

Although the velocity map of NGC 1052 and the central 
slit velocity curve superficially resemble a rotational veloc- 
ity field, the velocity structure shown in Figure [4] indicates 
that the NLR is actually dominated by an outflow. The color- 
coded velocity map shows only one of the two velocity com- 
ponents, where the component we plot shows redshifted and 
blueshifted sides that correspond to the jet-like features seen 
in the WFPC2 image. 

4.2. NGC 1497 

Due to the poor S/N of the spectra, especially from the 
first slit position, only a few rows in each exposure could be 
fit. Consequently, it is difficult to determine whether the gas 
is rotating. However, a velocity gradient of 150 km s" 1 is 
seen across 0"15 in the central slit position, and this could be 
consistent with rotation, although the velocity gradient is not 
symmetric about the location of the peak of the narrow-line 
region flux. 

4.3. NGC 2911 



Only a small number of velocity measurements could be 
made at each slit position due to the low S/N of the spectra. 
However, the velocity map suggests that the gas is dominated 
by rotational motion. The central slit position exhibits a gra- 
dient of 270 km s" 1 across the innermost 0"3. With higher 
quality spectra, this galaxy could be a good candidate for dy- 
namical modeling. The flux-weighted average velocity mea- 
sured in a 0"3 square aperture is 3285 km s -1 , while the sys- 
temic velocity taken from RC3 is 3167 km s -1 . Other optical 
measurements of the systemic velocity given by the NASA 
Extragalactic Database (NED) span a large range of veloci- 
ties, ~ 3030-3250 km s -1 . The inconsistencies between the 
various systemic velocity measurements may be the result of 
dust obscuration of the galaxy nucleus and one side of the 
emission-line disk, as a dust lane is visible in the WFPC2 im- 
age. 

4.4. NGC 3227 

The velocity field is very complex and chaotic. Two sepa- 
rate kinematic components are seen in all seven slit positions. 
At some locations, the velocity differences between the mul- 
tiple velocity components are as large as ~ 600 km s -1 . Our 
kinematic mapping of NGC 3227 shows that the NLR con- 
tains components with low velocities and low velocity disper- 
sions, as well as components with high velocities and high 
velocity dispersions. In particular, in the eastern direction in 
slit positions 1 and 2, one of the components has a velocity of 
~ 300 km s" 1 relative to the systemic velocity and a velocity 
dispersion of ~ 400 km s" 1 . In slit position 7, in the north- 
ern direction, one component has a velocity of ~ 300-350 
km s" 1 relative to the systemic velocity, and a velocity dis- 
persion of ~ 400 km s . Similar multicomponent structures 
have been seen in other Seyferts as wel l, including NGC 4151 
dKaiser et al.ll2000t iNelson et al.ll2000h . 

An additional Ha broad component was fit to spectra in all 
but the outermost two slit positions. The brightness of the 
Ha broad component in the central slit position, as well as 
in the adjacent slit positions 2 and 3, appears to be consistent 
with a bright point source that has been smeared by the STIS 
point-spread function (PSF). However, the flux of the broad 
component found at a Y Offset of -0"3 to 0."05 in slit position 
5, and the broad component flux at a Y Offset of -0"25 to 
0"15 in slit position 6, appear to be in excess of that expected 
from the wings of a point source. It is possible that NGC 3227 
exhibits scattered broad-line region light located in this region 
to the northwest of the nucleus. 

Previous studies of the ionized gas in NGC 3227 have 
also found chaotic velocity fields, as well as asym- 
metries and double-peaked profiles of the narrow Ha, 
[S II], [N II], and [O III] emission lines. These asym- 
metries were attributed to the superposition of sev- 
eral physically and kinematically distinct components, 
whe re the broader component was attributed to outflowing 
gas dArribas & M ediavillal ll994t iGonzalez D elgado & Perez! 
ll997HGarcia-Lorenzo et al.ll2001l) . Although strong outflows 
are seen in the ionized gas, the molecular gas in NGC 3227 
appears to be in regular rotation with a major axis at a position 
angle of 140° (ISchinnerereta l. 2000; Reunanen et alJl2002t 
IDavies et al.ll2006t iHicks & Malkanll2007l) . Our data show 
that one of the two velocity components, given by the black 
points in Figure [2d] and shown in Figure |U exhibits similar 
ov erall features as t h e mol ecular gas velocity field presented 
in IHicks & Malkanl (120071) . The redshifted and blueshifted 
ionized gas in this kinematic component roughly corresponds 
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spatially to the redshifted and blueshifted portions of the 
molecular gas. While one of the [N II] velocity components 
appears to trace the general features of the molecular gas, the 
ionized gas is still highly disturbed and is not participating in 
simple disk-like rotation. 

4.5. NGC4036 

The radial-velocity map shows the gas is rotating, although 
this motion is difficult to recognize from the multiple slit ve- 
locity curves alone. Velocities measured from the two slit 
positions to the southwest of the nucleus are blueshifted by 
as much as 200 km s -1 , while velocities from slit position 
2 are redshifted by 150 km s -1 relative to the systemic ve- 
locity of 1397 km s" 1 . The central slit position exhibits 
a velocity gradient o f ~ 300 km s" 1 across the inner 0."2. 
ICinzano et alJ (|1999[) examined the [O II] gas kinematics and 
found that NGC 4036 contained gas that was clearly not un- 
dergoing circular motion on scales less than 10". Addition- 
ally, based upon th eir WFPC2 narrow-band Ha+[N II] image, 
Pogge ^tail (l2000|) noted the existence of a complicated fila- 
mentary and clumpy structure. 

We did not require a Ha broad component in the fit to the 
spectra in any of the slit positions. How ever, through ground- 
based measurements. IHo et al.l d!997bl) determined that NGC 
4036 contained a weak, very low-amplitude, broad Ha com- 
ponent. This disagreement may be the result of variability, 
but another possibility is that the apparent broad Ha com- 
ponent in the ground-based data arises from the unresolved 
high-velocity material in the inner NLR. 

4.6. NGC 4278 

From the multislit velocity curves and the velocity map, 
the gas appears to be in organized motion and displays some 
characteristics of rotation. However, within the inner 0."5, 
the velocity field appears to be twisted, and may be the re- 
sult of a counter-rotating core. Another attractive possibility 
is that this central decoup led structure i s relat ed to the parsec- 
scale radio jet seen by iGiroletti et al.l (|2005). According to 
iGiroletti et al.1 d2005). the double-sided jet is seen at position 
angles of 155° and -40°, which suggests that the central kine- 
matic structure seen in the velocity field may be the result of 
material being pushed by the jet. The central slit position ex- 
hibits a velocity curve that is symmetric about the central row 
and has a gradient of ~ 400 km s" 1 across the inner 0."2. Two 
velocity components were needed to fit some of the spectra 
in the two outermost slit positions, where the additional com- 
ponent's velocity differed by ~ 450 km s -1 from the primary 
component. An additional broad Ha component was fit in the 
inner regions of the central slit position. 

4.7. NGC 4579 

Both the multislit velocity curves and the velocity map 
show that the gas is not in regular rotation. A few spectra in 
the off-nucleus slit position immediately adjacent to the center 
slit displayed two kinematically distinct components, where 
the velocity separation is ~ 450 km s . The nuclear disk 
may be nearly face-on since the primary velocity component 
is almost flat at ~ 1500 km s" 1 . In the central slit position, the 
[N II] velocity dispersion peaks at 440 km s -1 at a separation 
of 0"25 from the central row. A broad component of Ha was 
detected in the nuclear regions of the central slit as well as in 
slit position 4. In the STIS data, the nuclear spectrum reveals 
an asymmetric, very broad Ha component with a full width at 



zero intensity of ~ 18000 km s and "shoulders" on the red 
and blue side s dBarth et alj200l"ah . In their Ha+[N II] narrow- 
band image, Pogge et al.1 (120001) find complex, clumpy, and 
filamentary emission, which they liken to a shell or a ring. 
Ground-based measurements on larger scales also find irreg- 
ular kinematics (Go nzalez Delgado & Perezl ll996). The flux- 
weighted average velocity measured from a 1 '.'2 square aper- 
ture is 1465 km s , while the systemic velocity taken from 
RC3 is 1627 km s . The other optical measurements of the 
systemic velocity given by NED span a wide range of veloci- 
ties, -1480-1750 km s" 1 . 

4.8. NGC 4594 

The data for NGC 4594 show organized motion consistent 
with an overall rotation pattern, but with significant irregular- 
ities. The central slit position contains a velocity gradient of 
~ 550 km s" 1 across 0"25 Narr ow-band Ha+[N II] WFPC2 
images from Pogge et al. (2000) revealed a bright Ha core 
surrounded by spiral-like wisps th at extended to A" east and 
west and 1" south of the nucleus. Ems ellem & Ferruitl (120001) 
found that a cut through the [N II] A6583 velocity map along 
the galaxy's major axis showed the existence of a strong ve- 
locity gradient near the center, and that the gas kinematics 
within 1" was decoupled from the gas located in the spiral 
wisps. 

There has been some previous debate as to whether a broad 
Ha co mponent is present in NGC 4594. iKormendv et al.l 
( 1996) con cluded that such a comp onent existed; on the other 
hand, both lHo et al.l (Il997al) and lNicholson et all (1 19981) were 
able to reproduce the Ha+[N II] blend with only narrow com- 
ponents. In our analysis, an additional broad Ha component 
improved the fit to several of the nuclear rows of the central 
slit position. We also tried to fit multiple narrow components 
to each of the Ha, [N II], and [S II] lines without a broad com- 
ponent, but failed to find a satisfactory result. However, the 
Ha+[N II] lines are severely blended together and the S/N is 
low in the central rows, so we are unable to conclude that a 
broad component is definitely required. 

4.9. NGC 5077 

The velocity curves and the radial-velocity map show that 
the g as is likely to be in r egular rotation. This is in agreement 
with Ide Francesco et all (120081) . who used the same archival 
STIS dataset to perform gas dynamical modeling and deter- 
mined the black hole mass to be M B h = 6.8^g x 10 8 M . The 
flux-weighted average velocity measured from a square aper- 
ture of 0"3 is 2732 km s" 1 , while the systemic velocity taken 
from RC3 is 2832 km s" 1 . Other optical measurements given 
by NED cover the range - 2650 - 2850 km s" 1 . 

4.10. NGC 5635 

The velocity map and the velocity curve from the central 
slit position are consistent with disk rotation. The central slit 
position shows a velocity gradient of — 200 km s -1 across 
the nucleus. The spectra from the off-nucleus slit positions 
had such low S/N that only a single velocity measurement 
could be made. Unfortunately, the poor quality of the spectra 
precludes detailed kinematic modeling of the gas. 

4.11. NGC 6500 

The velocity field from the multiple slit positions is irregu- 
lar and not consistent with pure circular motion. Based upon 
the velocity map, the gas does appear to be dominated by rota- 
tion, though with some additional random components. This 
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is in agreement with lde Francesco et al.l (120081) . who also find 
complex velocity curves and departures from regularly rotat- 
ing gas, but they do note a general rotation trend from the 
central slit position. They suggest that the complicated ve- 
locity curves are the result of a nuclear expanding bubble. 
Gonzalez Delg ado"& Perezl (1996) detected multiple velocity 
components out to 1 kpc away from the nucleus, indicating 
that the extended NLR is experiencing an outflow. 

4.12. IC989 

The S/N is very low and velocities could only be mea- 
sured at a few positions, so the kinematic state of the gas is 
not well constrained. Wi th the same archival STIS dataset, 
Ide Francesco et all (|2008) note that the central slit position 
shows overall rotation with blueshifted gas southwest of the 
nucleus and redshifted gas northeast of the nucleus. They sug- 
gest the presence of a counter-rotating nuclear gas disk since 
the velocity gradient is reversed at the center. Alternatively, 
this may indicate the presence of a radial flow component. 

5. DISCUSSION 
5.1. Density Gradients 

In addition to determining the velocity, velocity dispersion, 
and flux from the measurement of emission lines, we mea- 
sured the [S II] A6716/[S II] A6731 and the [Ol] A6300/([S II] 
A6716+[S II] A6731) line intensity ratios. The [S II] ratio can 
be used to determine the electron density because of the dif- 
ference between the critical densities of the two lines, where 
Wen, = 1.5 x 10 3 cm" 3 for fS II] A6716 and n crit = 3.9 x 10 3 
cm" 3 for [S II] A6731 (e.g.. lPetersonlll997l) . However, this 
ratio is only sensitive to density over the range 10 2 cm" 3 s 
n e s 10 4 cm" 3 . The [O I]/[S II] ratio is also density-dependent, 
but not in a unique fashion since the ratio is influenced by 
the abundances, ionization level, and shape of the ionizing 
spectrum. However, the [O I]/[S II] ratio is sensitive to elec- 
tron density over a range wider than the [S II] ratio because 
of the larger critical densit y of the JO I] A6300 line, where 
«crit = 1-8 x 10 6 cm" 3 (e.g.. IPetersonll 19971) . Below, we show 
results from photoionization calculations to illustrate the use 
of the [O I]/[S II] ratio to probe high-density regions. 

Figure contains plots of the [S II] A6716/[S II] A6731 ra- 
tio, the [O I]/[S II] ratio, and the electron density determined 
from the [S II] ratio as a function of projected radial distance 
from the nucleus. The [O I] lines were not present in the NGC 
3998 spectra since the grating tilt was set to provide coverage 
of 6490-7050 A, so only the [S II] ratio and the electron den- 
sity are plotted as a function of distance from the nucleus. To 
examine the overall trends in electron density with radius, we 
computed average [S II] A6716/A6731 ratios in 0."l-wide ra- 
dial bins (or 0"05-wide bins for NGC 3998). With the radial 
averages of the [S II] A6716/[S II] A6731 ratio, the electron 
density as a function of radius was calculated using the IRAF 
task temden assuming a temperature of 10,000 K. In the in- 
stances where velocity splitting was present, the secondary 
component was too weak to produce an independent measure- 
ment of the electron density, and the two components were 
summed before the [S II] ratio was calculated. The remaining 
eight LINERs not included in Figure |5]have such large uncer- 
tainties in the density measurements that no information on 
the density structure could be extracted from the data. 

Significant density gradients are seen in NGC 1052, NGC 
4278, NGC 3998, and NGC 4579. Measurements of the den- 
sity structure in NGC 4579 have been previously presented 



by iBarth et al.l d2001al) . In these galaxies, the electron den- 
sity increases from ~ 200 - 900 cm" 3 in the outer regions to 
~ 1700- 16,000 cm" 3 at the center. In contrast to these four 
galaxies, NGC 4036 does not exhibit any systematic change 
in n e with respect to radial distance from the center or a cen- 
tral spike in density. It is interesting to note, however, that this 
galaxy did not need a broad Ha component to fit the central 
spectrum. Among the objects with a sufficient S/N to map n e 
from the [S II] lines, all of the broad-lined AGNs show central 
density peaks. 

The [O I]/[S II] ratio also rises dramatically toward the nu- 
cleus in some objects, increasing to as much as 1.2 ±0.1 in 
NGC 1052, 1.0 ±0.1 in NGC 4579, and 0.8 ±0.1 in NGC 
4278. These central [O I]/[S II] values are very large in com- 
parison to the average ratio typically found on scales of hun- 
dreds of parsecs from ground-based spectra of LINERs. For 
comparison, we used the emission-line measurements from 
iHo et al.1 d!997al) to calculate an average [O I]/[S II] ratio for 
both type 1 and type 2 LINERs. We first removed the galax- 
ies for which a "low confidence" or a "high uncertainty" flag 
was given for either the LINER classification or the [O I] /Ha 
and [S II] /Ha ratios. This left a sample of 57 galaxies from 
which we calculated the average [O I]/[S II] ratio to be 0.26 
in the ground-based 2"x4" aperture. From the same sample 
of 57 galaxies, those galaxies classified as type 1 LINERs had 
[O I]/[S II] ratios in the range 0. 15-0.54 with an average ratio 
of 0.32, while the type 2 LINERs had [O I]/[S II] ratios of 
0.13-0.58 with an average ratio of 0.21. 

Although the [O I]/[S II] ratios at the centers of NGC 1052, 
NGC 4278, and NGC 4579 are well above the typical val- 
ues found in LINERs, the [O I]/[S II] ratio is not a unique 
indicator of density. The ratio is additionally affected by the 
abundance ratio, the ionization parameter (U, which charac- 
terizes the number of ionizing photons per nucleon), and the 
shape of the ionizing continuum. In order to further interpret 
the high v alues of \Q I]/JS III , we used version 07.02.00 of 
CLOUDY dFerland et al.ll 998) to generate a few grids of sim- 
ple photoionization models. In particular, we used the AGN 
continuum model that takes as free parameters the tempera- 
ture of the thermal big blue bump component, the UV (2500 
A) to X-ray (2 keV) spectral slope a ox , and the power-law 
slopes of the low-energy part of the big blue bump (a uv ) and 
of the X-ray continuum (a x ). We left the thermal component 
temperature, a uv , and a x at their default values (T = 150,000 
K, a uv = -0.5, a x = -1). The a ox parameter was varied over 
-0.6 to -1 .6 in increments of 0.2, br acketing the typical range 
for low-luminosity AGNs (Ho 1999). Models were computed 
for a single plane-parallel slab of constant-density gas either 
at solar and twice solar metallicity, with gas density ranging 
from «h = 10 2 cm" 3 to 10 6 cm" 3 , and ionization parameter 
ranging from log U = -4 to -2. We calculated models that did 
and did not include dust grains. 

From these grids of photoionization models, we found 
small [O I]/[S II] ratios, consistent with ground-based results, 
at low values of «h (10 2 cm" 3 - 10 4 cm" 3 ). The ratio increases 
strongly for «h~10 4 cm" 3 . To illustrate our model results, in 
Figure|6]we plot the [O I]/[S II] ratio as a function of gas den- 
sity for the models with twice-solar metallicity and dust-free 
gas, and for a ox = -0.8 and -1 .4. For a ox = -0.8, a [O I]/[S II] 
ratio of 1 .0 corresponds to a density of ~ 2 x 10 4 cm" 3 for U 
= 10" 2 and 10" 4 , and a density of - 1 x 10 4 cm" 3 for U = 10" 3 . 
For a ox = -1.4, a [O I]/[S II] ratio of 1.0 corresponds to a den- 
sity of ~ 1 x 10 5 cm" 3 for U = 10~ 2 and 10~ 3 , and a density 
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of - 5 x 10 4 cm" 3 for U = 10" 4 . The smaller [O I]/[S II] ratio 
for a x = -1 .6 is expected since the increased UV flux relative 
to X-ray flux results in a less extended partially ionized zone. 
The models where solar metallicity was used and where dust 
was included showed qualitatively similar results to the mod- 
els presented in Figure [6] for the low ionization parameters, 
where densities greater than ~ 10 4 cm" 3 were needed in order 
to produce an [O I]/[S II] ratio of 1 .0. Thus, while the electron 
density cannot be directly determined from the [O I]/[S II] ra- 
tio, the rise in the line ratio seen in NGC 1052, NGC 4278, 
and NGC 4579 gives strong evidence for a steep inner rise in 
density. 

The single Seyfert galaxy in the sample, NGC 3227, ex- 
hibits a density gradient over a portion of the NLR, with 
n e ~ 600 cm" 3 in the outer regions and « e ~ 13, 000 cm" 3 at 
a radius of 10 pc. The electron density then decreases at the 
nucleus to ~ 2800 cm" 3 . Additionally, the [O I]/[S II] ratio 
rises toward the nucleus to 0.8 ± 0. 1 at a radius of 10 pc, then 
decreases to 0.4 ±0.1 at the center. 

Overall, the data from NGC 1052, NGC 3998, NGC 4278, 
and NGC 4579 provide clear evidence that Type 1 LINERs 
contain a radially density-stratified NLR. The detection of 
density gradients in Type 1 LINERs is consistent with pho- 
toionization models of LINERs. It also appears that density 
gradients can persist even in the presence of strong outflows 
such as those seen in NGC 1052 and those seen over a por- 
tion of the NLR in NGC 3227. Moreover, these five broad- 
lined AGNs exhibit similar density structure as can be seen 
in Figure [7] where we show the electron density as a func- 
tion of projected radial distance from the nucleus for the five 
AGNs. We fit a power law of the form n e = no(r/l pc)" to 
the density measurements from each galaxy separately. We 
calculated the average slope and normalization, as well as the 
standard deviation for each of the two parameters, and found 
«o = 5370 ±2495 cm" 3 and a = -0.60 ±0.1 3. Although there 
is a large range in the normalization values between the five 
broad-lined AGNs, the density measurements are character- 
ized by a similar slope. 

In order to determine whether unresolved density gradients 
exist even on the smallest angular scales probed by STIS, we 
examine the nuclear spectrum for evidence of a correlation 
between line width and critical density for collisional deex - 
citation (e.g., iFilippenko & Halpernlfl98lt lFilippenkol fT985). 
The line-width-n cr ; t relationship implies that the NLR has ra- 
dial gradients in velocity and density, where the broader emis- 
sion lines arise from gas closer to the black hole. In Figure[8] 
we plot the full width at half-maximum intensity (FWHM) of 
the [S II] A6731, [Nil] A6583, and [Ol] A6300 lines as a func- 
tion of critical density. The widths of these lines were mea- 
sured from the central row of the central slit position for all the 
galaxies in the sample with the exception of NGC 3227. Due 
to the sharp peak in the surface brightness at the nucleus of 
NGC 3227, the geometric rectification performed during the 
data reduction produced anomalous bumps in the continuum 
as well as an irregular [O I] A6300 line profile. We therefore 
binned together the innermost three rows from the central slit 
position before fitting Ha+[N II], the [S II] doublet, and the 
[O I] A6300 lines to a double-Gaussian narrow-component 
model. In the plot for NGC 3227 in Figure [8] we show the 
FWHM for each of the two Gaussian components for a single 
emission line. 

Of the eight objects in the sample with only [S II] and [N II] 
measurements, seven of the galaxies show that the FWHM of 
the [N II] line is greater than the FWHM of the [S II] line. 



Although many of these width measurements have large er- 
ror bars, it does appear that the FWHM rises with increasing 
critical density. Of the five galaxies with measurements of 
all three lines, NGC 1052 and NGC 4579 clearly show the 
presence of a line-width-n cr i t relationship, while NGC 4278 
exhibits an increase in the FWHM between the [S II] A6731 
and the [O I] A6300 line. NGC 4036 shows an increase in the 
FWHM between the [S II] A6731 and the [N II] A6583 lines, 
but a slightly smaller FWHM for the [O I] A6300 line. In 
contrast, NGC 3227 shows a decreasing FWHM with increas- 
ing critical density for both the high velocity dispersion and 
low velocity dispersion Gaussian components. The decreas- 
ing line width with increasing critical density indicates that 
the lowest density gas is actually moving the fastest, which 
may be the result of the strong outflow in NGC 3227. It is also 
possible, however, that this unusual trend seen in the Seyfert 
NGC 3227, as well as the decrease in the electron density seen 
in the inner region in Figure [5] may not be real and may be 
due to an inadequate fit to a rather complicated nuclear spec- 
trum. Without additional Seyferts in the sample, or higher 
luminosity AGNs similar to NGC 3227, it is impossible to 
determine whether this is an anomalous result or a more com- 
mon finding related to the properties of the AGN. With the 
exception of the Seyfert galaxy NGC 3227, a line-width-« cr ; t 
relationship can be observed even on the smallest STIS scales 
corresponding to ~ 8 pc, making the line-width-n cl i t relation- 
ship a generic property for LINERs. We are therefore seeing 
unresolved density structure, providing further evidence that 
the density continues to rise even within the innermost ~ 8 pc 
in most LINERs. 

5.2. Aperture Size and Line Width 

Another question that can be addressed with our measure- 
ments concerns the variations in emission-line width with 
aperture size. Such an analysis is useful in determining the 
scales over which the bulge dominates the NLR kinematics 
in low-luminosity AGNs. In the past, it has been assumed 
that ground-based observations sample NLR gas that is kine- 
mat ically dominated by the virial motion in the bulge p oten- 
tial ([N elson & Whittlel ll996l) . However, as discussed bv lLaorl 
(2003), low-luminosity AGNs have a more compact NLR 
whose kinematics may instead be dominated by the black 
hole. 

iRice et all (l2006h examined narrow-line widths as a func- 
tion of aperture size using single STIS slit observations of a 
sample of Seyfert galaxies. They found most of the galax- 
ies exhibited a constant or increasing FWHM with increas- 
ing aperture size. Our sample allows us to perform a similar 
analysis on low-luminosi ty AGNs. We a re able to consid- 
erably improve upon the IRice et all (120061) analysis by using 
the STIS multislit observations to study the line-width vari- 
ation within square apertures of increasing size, eliminating 
the biases introduced by studying only a relatively narrow and 
randomly oriented portion of the velocity field. 

For the six galaxies from program GO-7403, we choose 
square apertures of 0"2, 0"7, l."2, and l."7 on a side, and 
for the eight galaxies from program GO-7354 we used square 
apertures 0."1, 0."3, and 0."5 on a side. The sizes of the 
apertures were chosen such that the width of an entire STIS 
slit fell within the aperture. We then used the results of the 
spectral fitting to the narrow [N II] A6583 emission line and 
summed together the Gaussian components that fell within 
the various aperture bins. We characterized the profile of the 
summed Gaussians in each aperture by measuring the line 



9 



width at 50% of the peak height, thus defining a FWHM 
for the summed line profile. The FWHM was measured di- 
rectly from the summed model, not from a Gaussian fit to the 
summed model. 

The quantity FWHM/2.35 is plotted against half the length 
of the square aperture size (i.e., the "radius" of the aperture) in 
Figure [9] We plot FWHM/2.35 in order to compare with the 
bulge stellar velocity dispersion. Also included in the plots is 
an additional line-width measurement f or [N III A6583 f rom 
a ground-based 2" x4"aperture given bv lHo et all (Il997al) for 
all but three galaxies in our sample. The 2"x4"aperture is 
taken to be an effectively 3" square aperture and is plotted at 
a radius of l."5. Additionally, we found ground-based mea- 
surements of the bulge stellar velocity dispersion from the 
literature for 13 of the 14 galaxies in the sample. These 
values are listed in Table [T] and are plotted in Figure [9] as 
the dashed horizontal lines. We estimated the radius of the 
sphere of influence of the supermassive black hole, given by 
r g = GMbh/o",, in each of the 13 galaxies for which a stellar 
velocity dispersion m easurement could be found. We used the 
iTremaine et al.l d2002l) Mbh-c* relation to estimate the black 
hole mass in all of the galaxies except NGC 3227, NGC 3245, 
NGC 3998, NGC 4594, and NGC 5077. For these galaxies we 
used black hole mass measurements from the literature (Table 
[T). The sphere of influence was then converted into an angu- 
lar scale using the distances listed in Table Q] and plotted as 
the dotted vertical lines in Figure [9] 

For most of the galaxies in the sample, the [N II] veloc- 
ity dispersion peaks within the black hole sphere of influence 
and then decreases in larger apertures, approaching the bulge 
stellar velocity dispersion. This implies that on scales acces- 
sible to ground-based measurements, the width of the low- 
ionization lines, such as [N II], is set by the bulge and not 
by the black hole. The work by lLaorl ([2003) suggested that 
in low-luminosity AGNs, the NLR may be dominated by the 
black h ole and not the bulge. However, the arguments put 
forth by |Laorl(l2003l) involved the [O III] line, which due to the 
higher ionization and higher critical density, originates from 
a region closer to the black hole than the [N II] line. Also, 
our finding that, for most galaxies in the sample, the veloc- 
ity dispersion is largest within the black hole sphere of influ- 
ence supports the idea that STIS central line-width measure- 
ments trace the black hole mass, and can therefore be used 
to place upper limits on black hole m asses in nearby galaxies 
dSarzi et al.ll2002t iBeifiori et al.ll2007l) . 

NGC 1052, NGC 3227, and NGC 3245, however, show 
a different trend, where the FWHM increases with aperture 
size. Both NGC 1052 and NGC 3227 contain successively 
wider or roughly constant [N II] narrow lines for the increas- 
ing STIS aperture sizes, but a smaller FWHM for the ground- 
based aperture. In the cases of NGC 1052 and NGC 3227, the 
strong radial motions create large velocity splitting between 
components, causing an increase in the FWHM. The smaller 
line width seen in the ground-based aperture may be the result 
of the outflows decelerati ng at larger radii, as describe d in the 
outflow cone models bv lCrenshaw & Kraemerld2000l) . 

An initial expectation for the behavior of line width with 
aperture size for gravitationally dominated gas may be a de- 
creasing FWHM with increasing aperture size. This is reason- 
able since most of the ionized gas flux comes from a region 
close to the black hole, where the gas is moving the fastest 
and thus has the largest line width due to unresolved rotation. 
However, the line-width variation with aperture size in NGC 
3245 clearly disagrees with the above expectation, since the 



gas in NGC 3245 is know n to be in near-perfect disk-like ro- 
tation (Barth et al. 200Tbj). In order to gain a more complete 
understanding of line-width variation with aperture size, we 
constructed a toy thin-disk model for gas in circular rotation. 
First, we generated a model velocity field on a grid where 
each model pixel was 0"00507 on a side (i.e., 10 x oversam- 
pled relative to the STIS pixel size). The velocity field is set 
by the black hole mass, the disk inclination a ngle, and the 
stellar mass profile. We used a iHernquistl (1 1 9901) model to de- 
scribe the stellar mass contribution. The intrinsic line-of-sight 
velocity profiles were assumed to be Gaussian before passing 
through the telescope optics, and were calculated on a veloc- 
ity grid with a bin size of 25.24 km s" 1 , in order to match the 
STIS pixel scale. 

The model velocity field was then projected onto the plane 
of the sky, and was synthetically "observed" in a manner that 
closely matches the STIS observation. This synthetic observa- 
tion involved summing the intrinsic line-of-sight velocity pro- 
files at every point in the disk and including the blurring intro- 
duced by the nonzero STIS pixel size and the telescope PSF 
In particular, the line-of-sight velocity profiles were weighted 
by the emission-line surface brightness profile, modeled by 
the function S(r) = So + S\e~ r / r °. We included an intrinsic ve- 
locity disp ersion, similar to those observed in oth er early-type 
galaxies dvan der Marel & van den Bosch 1998), by setting 
the widths of the line-of-sight velocity profiles according to 
<r(r) = <7o + o\e~ r l r " . Fina lly, we accounted fo r the telescope 
PSF by using Tiny Tim dKrist & Hookll2004T) to generate a 
0"2x0."2 subsection of the full STIS PSF. The PSF was sub- 
sampled to match the model pixel grid and computed for a 
monochromatic filter passband at 6600 A. Each velocity slice 
of the pixel grid was convolved with the PSF. These model- 
ing steps ar e a simplified versi on of the modeling described 
in detail by iBarth et al.l (|2001bh . For our purposes, this syn- 
thetic observation involved placing a single STIS slit with in- 
creasing aperture lengths at different angles mea sured from 
the ga seous disk major axis in order to mimic the iRice et all 
(2006) analysis. We additionally placed square apertures of 
increasing size over the model velocity field to replicate our 
measurements previously discussed. In both instances, we 
summed the line profiles that fell within the various apertures 
and calculated the FWHM of the summed profile. 

As an example, in Figure [10] we show the results of the 
toy model using the properties of NGC 3245 and NGC 3998 
as the input parameters. We selected these galaxies because 
in both cases the gas is known to be in circular rotation but 
the behavior of line width versus aperture size is quite dif- 
ferent for the two objects. We used the black ho le mass and 
inclinat ion a ngle determined from the models of Bart h et all 
d2001bl) a n d Ide Francesco et a l. (2006), and we choose the 
iHernquisIl (1990) bulge normalization and the scale radius 
such that the toy model velocity field would match the ob- 
served velocity field as closely as possible for the outermost 
velocity measurements. Similarly, we determined the param- 
eters of the emission-line surface brightness profile by requir- 
ing the profile to match the observed [N II] flux distribution. 
Also, for NGC 3245, the parameters of the function describ- 
in g the widths of the l ine-of-sight velocity profiles were given 
bv lBarthetai](l2001bl). 

Note that Ide F rancesco et al.l d2006l) did not include any 
intrinsic velocity dispersion in their model of NGC 3998. 
We therefore first generated disk models using only a small 
constant intrinsic line width of ctq = 20 km s" 1 . We found 
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an overall trend that qualitatively matched the data, but the 
line widths overall were too small, as can be seen in Fig- 
ure [TO] and the line profile was double peaked in each of 
the apertures. Since a double-peaked line profile is not seen 
in the data, we instead included an intrinsic velocity dis- 
persion with the parameters <tq = 220 km s" 1 , <j\ = 100 km 
s -1 , and r„ = 10 pc. We chose this intrinsic velocity disper- 
sion since it produced a central line width more consistent 
with the data and a single-peaked line profile in each of the 
apertures. Such values for the intrinsic velocity dispersion 
are not unreasonable given previous r esults of dis k modeling 
by Ivan der Marel & van den Boschl (119981) and iBarth et all 
(2001b). Moreover, we stress that we are employing a sim- 
plified version of the disk modeling; a complete treatment, in- 
cluding a more appropriate emission-line surface brightness 
profile and a more realistic model of how the light passes 
through the telescope optics, would change the values of the 
intrinsic velocity dispersion parameters needed to reproduce 
the observed line widths. 

The top panels for both NGC 3245 and NGC 3998 in Figure 
[TOlare the results from using progressively larger square aper- 
tures. Since we adopted a simplified disk model, the models in 
Figure QJJdo not exactly numerically match the observational 
data. However, with our model, we are able to reproduce the 
qualitative behavior of line-width variation with aperture size. 
The model shows that the line width increases with increasing 
aperture size for NGC 3245, and the FWHM decreases with 
increasing aperture size for NGC 3998. The additional plots 
in each column are the results of using a single STIS slit, at 
an orientation that is aligned along the major or minor axis of 
the gaseous disk. Clearly, the orientation of the single slit not 
only affects the numerical values of the line widths, but also 
the general behavior of the line-width variation with aperture 
size, as can be seen in NGC 3245. In both cases, the disk 
models always produced a monotonic trend with line width, 
either an increasing or decreasing line width with increasing 
aperture size. 

These models demonstrate how either an increasing or de- 
creasing line width with increasing aperture size can be seen 
for gravitationally dominated gas. NGC 3245 contains a mas- 
sive bulge that produces large circular velocities at distances 
far from the nucleus. While this is also true for NGC 3998, 
NGC 3245 has a flatter emission-line surface brightness pro- 
file with relatively more flux at larger radii than in NGC 3998, 
as can be seen in Figure [3] Additionally, the gaseous disk in 
NGC 3245 is highly inclined at i = 63°, and the disk in NGC 
3998 is less inclined with i = 30°. The combination of the 
flatter surface brightness profile, the highly inclined disk, and 
the massive bulge of NGC 3245 results in the contribution 
of a significant amount of gas with large circular velocities 
further away from the nucleus. Consequently, the line width 
increases with increasing aperture size. 

Therefore, based upon the results of the toy model and the 
analysis of the observational data, the overall behavior of line 
width with aperture size is not a unique indicator of the kine- 
matic state of the gas. Outflowing gas can give an increas- 
ing line width with aperture size, as seen in NGC 1052 and 
NGC 3227, but decelerating outflows could produce a de- 
cre asing trend in lin e width with aperture size as discussed 
by iRice et alj (120061) . Similarly, rotation can generate either 
a decreasing or increa sing line w i dth wi th increasing aperture 
size. Thus, while the IRice et aD (|2006) results are useful for 
identifying which aperture sizes yield similar stellar and gas 
line widths, our results show that measurements such as these 



are not necessarily useful as an indicator of whether the line 
widths are primarily virial in origin. 

Finally, many previous studies have examined the possible 
role of outflows and nongravitational motions to explain the 
scatter in the <J g -cr* relationship. However, our modeling and 
data have shown that, even for galaxies whose NLR gas is in 
overall rotation and dominated by virial motion, the gas line 
widths can differ significantly from the stellar velocity disper- 
sion, as seen in NGC 3245 and NGC 3998. We find that the 
gaseous disk inclination, the surface brightness distribution, 
and the aperture size and shape are important factors that also 
most likely contribute to the scatter seen in the <7g-er* relation- 
ship. 

6. CONCLUSIONS 

We have used STIS data to measure detailed velocity, veloc- 
ity dispersion, flux, and line-intensity ratio information within 
~ 100 pc from the nucleus of 14 nearby galaxies. All of the 
galaxies harbor low-luminosity AGNs, 13 of which are classi- 
fied as LINERs and one of which is a Seyfert. We see a large 
diversity in the velocity fields of the 14 galaxies, where simple 
disk rotation is not a common state. Only a few galaxies have 
definite disk rotation, while some others have possible disk ro- 
tation but the emission-line S/N is too low to clearly determine 
the kinematic state of the gas. Some velocity fields are domi- 
nated by irregular motions, some show an overall rotation but 
with large random components, and others clearly exhibit out- 
flows. Dynamical modeling of NGC 291 1 and NGC 4594 to 
determine black hole masses may be possible and worthwhile, 
although chaotic and disorganized motion will limit the accu- 
racy of mass measurements. 

Through the [S II] A6716/[S II] A6731 diagnostic, we were 
able to detect significant electron-density gradients in four 
(NGC 1052, NGC 3998, NGC 4278, and NGC 4579) of the 
five LINERs having a sufficient S/N. Moreover, the electron- 
density gradients in these objects share a similar shape that 
can be described by a power law with a slope of -0.60 ±0.1 3. 
A similar increase near the nucleus was seen in the [O I]/[S II] 
line ratio for the first three of these galaxies. These results 
demonstrate that density gradients within the inner ~ 20 pc 
are common in type 1 LINERs. This is consistent with ex- 
pectations from photoionization models of LINERs and from 
ground-based detections of correlations between line width 
and n cr j t in LINERs. We find that the density gradients persist 
in NGC 1052 and over a portion of the NLR in NGC 3227, 
even in the presence of strong outflows that clearly have a 
large impact on the NLR gas as seen through the disorganized 
velocity fields with multiple velocity components. We also 
see evidence for correlations between line width and n cl ; t in 
most LINERs even on the 0"1 scales probed by STIS. 

Measurements of the STIS spectra additionally provided 
a means to study the [N II] A6583 narrow-line width varia- 
tion with spectrograph aperture size. For most galaxies in our 
sample, the emission-line velocity dispersion peaks within the 
black hole sphere of influence, then decreases in larger aper- 
tures, approaching the bulge stellar velocity dispersion. This 
trend is consistent with, but does not uniquely imply, virial 
motion. The increasing line width with aperture size seen in 
NGC 1052 and NGC 3227 is due to the outflows that domi- 
nate the NLR kinematics. In NGC 3245, the increase is the re- 
sult of the combination of the galaxy's fairly flat emission-line 
surface brightness profile, the galaxy's massive bulge, and the 
highly inclined gaseous disk. NGC 3245 demonstrates that 
even gas in pure disk rotation can exhibit an increasing line 
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width as a functio n of aperture size . 

While work by iBorosonl (120051) and iGreene & Hoi (120051) 
provide evidence for a connection between the Eddington ra- 
tio and kinematic disturbance in the NLR, which may reflect 
an underlying trend of outflow-dominated NLRs being more 
predominant at high AGN luminosities, our sample of galax- 
ies is too small and does not span a sufficiently wide lumi- 
nosity range to arrive at such a conclusion. Since such a re- 
lationship would be helpful in understanding AGN feedback 
and the effects on host-galaxy properties, further investiga- 
tion, ideally with IFU data on a much larger sample of AGNs 
spanning a wide range of luminosities, is needed. 
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TABLE 1 

Properties of Galaxies in the Sample 



Galaxy 


Morphological 


cz 


D 


Spectroscopic 






log L(Hq) 


Referenc 




Type 


(kms- 1 ) 


(Mpc) 


Classification 


(km s-') 


(M@) 


(ergs s~' ) 


for Oi, 


NGC 1052 


E 


1474 


14.7 


LI. 9 


215 


1.8 x 10 


39.28 


1 


NGC 1497 


SO 


6170 


80.8 


L2* 


249 


3.3 X 10 8 


39.05 1 " 


2 


NGC 2911 


SO pec 


3167 


45.1 


L2 


237 


2.7 X 10 8 


39.44 


3 


NGC 3227 


SAB pec 


1145 


18.3 


SI. 5 


128 


1.5 x 10 7a 


40.26t 


1 


NGC 3245 


SO 


1348 


20.3 


T2: 


232 


2.0 x 10 8b 


39.51 


4 


NGC 3998 


SO 


1066 


16.4 


L1.9 


333 


2.7 x 10 8c 


39.76 


1 


NGC 4036 


so 


1397 


21.5 


L1.9 


166 


6.7 x 10 7 


39.23 


5 


NGC 4278 


E 


628 


13.1 


L1.9 


270 


4.5 x 10 8 


39.43 


6 


NGC 4579 


SABb 


1627 


16.9 


S1.9/L1.9 


160 


5.5 x 10 7 


39.44 


1 


NGC 4594 


Sa 


1082 


19.4 


L2 


257 


1.0 x 10 9d 


39.67 


7 


NGC 5077 


E 


2832 


AAA 


L1.9 


275 


6.8 x 10 8e 


39.75 


8 


NGC 5635 


S pec 


4352 


64.3 


L2* 






38.91 1 " 




NGC 6500 


Sab 


2975 


44.6 


L2 


168 


6.7 x 10 7 


40.42 


9 


IC989 


E 


7556 


110.2 


L2* 


176 


8.1 x 10 7 


39.47t 


10 



References. — (1) Nelson & Whittle 1995; (2) Wegner et al. 2003; (3) di Nella et al. 1995; (4) Tonry & Davis 1981: (5) Fisher 1997; 
(6) |BenderetaU1994t (7) 1 van der Marel et all 1991 (8) |Faber et all 19891 (9) |Terlevich et all 19901 Q0) IWhiteet all 19831 

NOTE. — Morphological types and distances are taken from the NASA Extragalactic Database (NED), where the distances are derived 
using the Mould et al. 1 200G) Virgo + Great Attractor + Shapley Supercluster Infall velocity field model. Recession velo cities come from 
the optical line measure ment of the Third Reference Catalogue of Bright Galaxies (RC3) (He Vaucouleurs et al. 1991). Spectroscopic 
classifications come from Ho et al. iT997jj), where S = Seyfert nucleus, L = LINER, T = Transition object, and the number next the class 
specifies the type. Type 2 objects have no detectable broad lines, type 1 .9 objects show a broad Ha component, and type 1 .5 objects exhibit 
broad Ha and H/3 lines. An uncertain spectroscopic classification is given by The spectral classifications mar ked with an asteri sk were 
based upon whether a broad Ha component was seen i n the data from program GO-7354 and were not taken from[Ho et al. 1 1997$. Black 
hole masses, unless marked, were estimated using the Tremaine et al. 1 2002.) Mbh — c* relation. The black hole masses for NGC 3227, 
NGC 3245, NGC 3998, NGC 4594, and NGC 5077 were determined through dynamical measurements. Ha luminosities were calculated 
using the distances listed in column (4) and the Ha flux taken from, Ho et al. 1 1997a) and Ho et al. ( 2003), except for those objects marked 
with a cross: NGC 1497, NGC 3227, NGC 5635, and IC 989. The Ha luminosities for these four objects were calculated from the flux 
measured within the largest aperture that could be constructed given the number of STIS slit positions for the galaxy and the number of rows 
from which spectra could be extracted. 

a |Davies et ail 2006 . b Barth et al. 2001b. c fdeFrancesco et alT2006. cl Kormendy et al. 1996. e de Francesco et al|2008l 
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Summary of Observations 
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of 
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(s) 
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Note. — 



The exposure time column gives the average of each slit position's exposure time for each galaxy. 
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FIG. 1. — Examples of various models used to fit the spectra: (a) a spectrum from NGC 4036 with a single-Gaussian fit to the narrow lines, (b) a spectrum 
from NGC 3227 having multiple velocity components with double-peaked profiles and a broad Ho component, and (c) a central spectrum from NGC 3998 fit 
with single Gaussians to the narrow lines and a broad Ha component. Thin solid lines represent the individual Gaussian components, and the sum of the model 
components is given by the thick solid line. 
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FIG. 2a. — The velocity, velocity dispersion, and flux as a function of position along the slit for each of the slit positions in NGC 1052. Spectral rows that were 
binned together in order to improve the S/N are shown with error bars in the ^-direction, where the error bars denote the radial range over which the spectra were 
extracted. The grey points represent separate velocity components. The horizontal dashed line shows the galaxy's systemic velocity that was taken from optical 
line measurements in RC3. The positive and negative y-offsets correspond to the spectra extracted from the rows above and below the central row, respectively, 
when the slits are oriented such that the STIS instrumental y-axis points directly upward. Thus, the bottom half of the slit drawn over the F658N WFPC2 image 
corresponds to the left side of the velocity, velocity dispersion, and flux plots. The length of the slits in the same diagram corresponds to the length over which 
we could measure the emission lines. The upper-right panel shows the spectrum extracted from the nuclear region of the central slit, where the flux density is in 
units of lCT 18 ergs s _1 cirr 2 A -1 . No WFPC2 continuum images were available. 
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FIG. 2b. — NGC 1497; see Figure l2alfor description. 
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FIG. 2c— NGC 291 1; see Figure l2alfor description. 
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FIG. 2d.— NGC 3227; see Figure|2a]for description. 
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FIG. 2e. — NGC 3245; see Figure[2a]for description. The S/N was too low to accurately fit the [S II] lines. See also lBarth et alJKOOlbl) . 
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FIG. 2f. — NGC 3998; see Figure[2a]for description. See also de Francesco et al] I200fl) . 
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FIG. 2g. — NGC 4036; see Figure|2a]for description. 
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FIG. 2l. — NGC 4579; see Figure[2a]for description. The spectra near the central regions in the second slit position contained severely blended Ha and [N II] 
lines, where the widths of the three Gaussians were held fixed at a specific value in order to achieve the most reasonable fit. The open circles correspond to rows 
where this width restriction was necessary. 
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FIG. 2j. — NGC 4594; see Figurel2alfor description. 
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FIG. 2k. — NGC 5077; see Figure|2a]for description. See also de Francesco et al. 12008). 
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FIG. 2l. — NGC 5635; see Figure [2a]for description. 
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FIG. 2m.— NGC 6500; see Figure[2a|for description. See also de Frances coeTaTl HOPS) . 
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FIG. 2n. — IC 989; see Figurel2alfor description. See also de Francesco et al. (2008). 



29 



m 
u 

0) 



CO 
CO 
LO 
CD 

< 



900 
600 
300 




0.5 



I 1 1 1 1 I 
NGC 1052i 



1 1.5 



FT 



100 
50 




NGC 291 1 - 



0.1 



0.2 



100 -' 



50 




I 1 1 1 1 I 1 L 
NGC 3245 



"I i i i i I 
0.5 



rr 

_m 

800 1 
400 -> 



NGC 4579- 



1 1.5 




NGC 3998 



-I ^for*j » , ! ■ | l ») *l| M ^ l h'i -1 







2000 
1000 



0.5 



1.5 



' 1 1 I 1 1 - 
NGC 4594: 



5* 



200 FT 1- 
- _s 

150 E- 
100 '- m M 

50 L 
Bl 



0.2 0.4 0.6 

' i I i ' i ' I ' i. 

NGC 4036: 



|-i i i . f*?*f A fT1*f*i i * 

0.2 0.4 0.6 



T 



0.5 



600 
400 r 
200 r 

o El 
o 



NGC 50771 



0.1 0.2 



0.3 




_ 1 1 1 1 


i i i i | i i i _ 










i i i | i 




NGC 4278: 


200 








NGC 56351 


- I ^ m m 














- ^jaB 




100 




t 

i 


r 

L 




- ME H^*"K 














11 , , , 









i i 




1 , , , , f , i 



400 



200 - I 



I 1 1 1 I 1 
NGC 6500 



0.5 1 

120 FT 
80 7T 
40 '-. 

o Bl 



0.1 



0.2 



I 1 1 1 1 I 1 ±| 
IC 989 



0.2 0.4 0.6 0.1 0.2 0.3 

Projected Distance (Arcsec) 

FIG. 3. — [N II] A6583 surface brightness as a function of projected angular distance from the nucleus. In the plots of the galaxies from program GO-7403, 
the black triangles are averages of the measurements in 0" 1 bins, and the grey squares represent the entire, unbinned set of measurements. In the plots of NGC 
3998, NGC 4594, and NGC 6500, the black triangles are averages of the surface brightness within a 0"05 bin. 
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FIG. 4. — Smoothed radial-velocity fields in the slit frame generated by combining the multisl it veloci ty measurements for each galaxy. For the galaxies having 
multiple velocity components, we show the velocity component plotted in black in Fi gures l2al 1 2dl l2hl and|2i] The minimum and maximum velocities for each 
galaxy relative to the systemic velocity are given next to the key. For all of the galaxies, green corresponds to the systemic velocity. Since the galaxies from 
program GO-7403 were observed with slits separated by 0f'25, we interpolated over the spacing to determine the velocities at those locations. 
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FIG. 5.— [S II] A6716/[S II] A6731 line ratio, [O I] A63()0/([S II] A6716+[S II] A6731) line ratio, and electron density as determined from the [S II] line 
ratio, as a function of projected radial distance from the nucleus. The black triangles are averages of the measurements in 0"1 bins, and grey squares represent 
the entire, unbinned set of measurements. In NGC 3998, the black triangles are averages of the ratios within a 0"05 bin. 
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FIG. 6. — [O I] A6300/([S II] A6716+[S II] A6731) as a function of density for a range of values of the ionization parameter. These photoionization models 
were generated using CLOUDY with a metallicity of two times solar in the absence of dust and with the AGN continuum model where ct ox = —0.8 and ct ox = —1.4. 
The models where solar metallicity was used and where dust was included show qualitatively similar results for the low ionization parameters. 
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FIG. 7. — Electron density determined from the [S II] line ratio as a function of projected radial distance from the nucleus for the five galaxies exhibiting 
a significant density gradient. Error bars have been omitted for clarity, but are given in Figure \5\ NGC 1052 (blue filled squares), NGC 3227 (purple filled 
triangles), NGC 3998 (orange crosses), NGC 4278 (red filled circles), and NGC 4579 (green open squares) show similar electron-density gradients, where the 
best-fit power law, with a slope of —0.60 ± 0. 13, is plotted as the dashed black line. 
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FIG. 8. — FWHM of the [S II] A6731, [N II] A6583, and [O I] A6300 lines measured from the central row of the central slit position for each of the galaxies as 
a function of the critical density, with the exception of NGC 3227 where the innermost three rows were binned together. The correlation between line width and 
n cri , is seen in most of the LINERs in the sample, even on the Of 1 scales probed by STIS. For NGC 1052, NGC 3227, NGC 4278, and NGC 4579, the error bars 
are smaller than the plot symbols. Both the low velocity dispersion (plotted in black) and high velocity dispersion (plotted in grey) components for NGC 3227 
are shown. 



34 




FIG. 9. — [N II] A6583 line width as a function of aperture size. Depending on how many slit posit ions were used t o observe the galaxy, either two, three, 
or four aperture bins were constructed. The outermost point comes from a ground-based measurement 1H0 et al. 1997a) using effectively a 3" square aperture. 
NGC 1497, NGC 5635, and IC 989 were not included in the Ho et al. ( 1997a) sample and thus no measurement for l'.'5 is given. The abscissa represents half 
of the square aperture size in arcseconds. The dotted vertical line is the black hole sphere of influence and the dashed horizontal line is the bulge stellar velocity 
dispersion measured from a ground-based aperture. 
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FIG. 10. — FWHM/2.35 as a function of aperture size as determined through the use of a circularly rotating thin-disk toy model for NGC 3245 and NGC 3998. 
The black hole sphere of influence for each galaxy is plotted in the dotted vertical line and the bulge stellar velocity dispersion for each galaxy is shown by the 
dashed horizontal line. The variation in line width with aperture size in the top set of plots was found using successively larger square apertures. For NGC 3998, 
the black filled squares are the result of using a disk model with a radially dependent intrinsic velocity dispersion, while the grey open squares are the result of 
using a model with a small constant intrinsic velocity dispersion. The bottom two sets of plots were generated by sampling the modeled velocity field with a 
single slit of width Q'.'2 for NGC 3245 and Q'.'l for NGC 3998, with progressively longer lengths. The single STIS slit was aligned along the major axis of the 
gaseous disk for the middle set of plots, and along the minor axis for the lower panels. The abscissa represents half of the aperture size in arcseconds. 



